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Abstract 



We have received rigorous relations between the exact and first- 
order parameters of the Moliere multiple scattering theory, instead of 
the approximate one obtained in the original paper by Moliere. We 
also estimated the corresponding Coulomb corrections to the Born 
approximation of the parameters and studied their Z dependence in 
the range of nuclear charge from Z = 4 to Z = 82. Additionally, 
we evaluated absolute and relative accuracies of the Moliere theory in 
determining the screening angle and have concluded that for Z ~ 80 
they are about 20 and 34 percents, respectively. 
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1 Introduction 

The Moliere multiple scattering theory of charged particles' [TJ [2] is the most 
used tool for the multiple scattering effects accounting at experimental data 
processing. The experiment DIRAC [5] and many others [4-7] meets the 
problem of the excluding the multiple scattering effects in matter from ob- 
tained data. The standard theory of multiple scattering proposed by Moliere 
and some its modifications [8-11] are used for this aim. 

Estimation of the theory accuracy is of especial importance in the case 
of DIRAC experiment for it's high angular resolution. One of the significant 
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sources of the inaccuracy of the Moliere theory is the usage in [H [2] an 
approximate relation for the exact and first-order values of the screening 
angle (% a ) obtained in the original paper of Moliere pQ: 

Xa ~xf^l + 3.34(Za//3) 2 , (1) 

since this parameter also determined other important quantities of the theory. 
In this connection, the problem of receiving the corresponding exact relation 
and estimating the Moliere theory accuracy becomes important. 

In the present work, we have obtained rigorous relations for a number 
of exact results of the theory and their Born approximations. We have also 
evaluated numerically the Coulomb corrections to the Born results and stud- 
ied the Z-dependence of these corrections. In addition, we have estimated 
the absolute and relative accuracy of the Moliere theory in determining the 
screening angle Xa- 

The outline of the paper is as follows. In Sec. 2-4 we review some basic 
results of [H |2]. Namely, a solution of the transport equation (Sec. 2), the 
Moliere expansion for this solution (Sec. 3), and determining the screening 
parameters by Moliere (Sec. 4). The results of the present work are obtained 
in Sec. 5-6. In Sec. 5 we consider an another determination of the screen- 
ing parameters allowing to receive rigorous relations between their exact and 
first-order values. In Sec. 6 we calculate the numerical values of the Coulomb 
corrections to the Born approximation in the range Z = 4 to Z = 82 and es- 
timate the inaccuracy of the Moliere theory. Finally, in Sec. 7 we summarize 
the main results of the paper. In Appendix we present an alternative way of 
obtaining the approximate solution of the transport equation for the thick 
targets. The some preliminary results of this work were announced in Refs. 



2 The transport equation and its solution 

Let all scattering angles are small 9 <C 1 so that sin 9 ~ 9, and the scattering 
problem is equivalent to diffusion in the plane of 9. Now let cr(x) be the 
elastic differential cross section for the single scattering into the angular 
interval x = 9 ~ ®' i an d W M (9,t)9d9 is the number of scattered particles 
in the interval d9 after traversing a target thickness t. Then the transport 
equation is 

9WM ^ t] = -n W u (9,t) I a{ X )d 2 x + n I W M {9 - x,t)o{ X )d\, (2) 
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where uq = cNa/A is the number of the scattering atoms per cm 3 , Na is the 
Avogadro number, and A is the atomic weight of the target atoms. 

Moliere solved this equation for the determination of the spatial-angle 
distribution function W M (9,t) by the Fourier-Bessel method and gets to her 
a general expression 



W M (0,t) = J J (67j)g{ri,t) V dn, (3) 

Q 

in which 

g( V ,t) = exp[N(r ] ,t)-N (0,t)}, (4) 

9 is the polar angle of the track of a scattering particle, measured with respect 
to the initial direction z\ rj is the Fourier transform variable corresponding 
to 6; and the Bessel function J of order is an approximate form for the 
Legendre polynomial appropriate to small scattering angles f2\ IT5]. 
In the notation of Moliere 



N(r], t) = 2im t j cr(x)Mxv)xdx, (5) 
o 

and iVo is the value of ^ for rj = 0, i.e. the total number of collisions 

oo 

N (0,t) = 2nn t [ a( X )xdx- (6) 



The magnitude of No — N is much smaller than Nq for values rj, which are 
important, and may be called the 'effective number of collisions'. 
Inserting Eqs. (|4])-(|6]) back into (J3j) , we have 



W M (6,t) = / rjdr)J Q (9rj)exjp 



-2nn t / ct(x)xg?x[1 - J (xv)) 



(7) 



This equation is exact for any scattering law, provided only the angles are 
small compared with a radian. 

At g(r), 0) = 1 for all rj the expressions (j3])-([6]) can be rewritten as follows: 

oo 

W M (9,t) = [ J (9rj)e- n ° tQ ^rjdrj, (8) 
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where 



g(r 7 )=27r / <r(x)[l-Mxv)]xdx- (9) 



o 



This result is mathematically identical with result of Snyder and Scott for 
the distribution of projected angles [14] . 

3 The Moliere expansion for the transport 
equations solution 

One of the most important results of the Moliere's theory is that the scat- 
tering is described by a single parameter, the so-called 'screening angle' (x a 

Of X'aY 

X' a = VU07 X a = [exp (C E - 0.5)] X a ~ 1-080 X a, (10) 

where C E = 0.57721 is the Euler's constant. 

More precisely, the angular distribution W M (8)6d6 depends only on the 
ratio of the 'characteristic angle' Xc which describes the foil thickness, to the 
'screening angle' which describes the scattering atom: 

b = \n(^] =ln(^) + l-2C B ~lnJV . (11) 



Xa / V Xa , 

The screening angle Xa can be determined approximately by the relation 
2 ~Xo(l.l3 + 3.76 a 2 ) = ( X f) 2 (l + 3.34a 2 ) (12) 



X. 

with the so-called 'Born parameter' 

Ze 2 Za 



hv (3 



(13) 



The second term in (12) represents the deviation from the Born approxi- 
mation. If the value of this term is 0, the value of the screening angle is 

Xa = Xa = XO VTI3. 

The angle xo is defined by 

Z l / 3 m e Z l / 3 m e a , s 

y = 1.13 -= — , 14 

AU 137p 0.885p ' V ; 

where p = m e v is momentum of incident particle and v is its velocity in the 
laboratory frame. 
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The characteristic angle is defined as 

Jp) ' (15) 

Its physical meaning is that the total probability of single scattering through 
an angle greater than Xc is exactly one. 

Putting XcV = V and setting 9/xc = u, for the most important values of 
rj of order l/xc, we get Moliere's transformed equation 



W M (9)9d9 = udu I ydyJ (uy)exp< — 



b — In 



(16) 



which is very much simpler in form than ([7]). 

In order to obtain a result valid for large all angles, Moliere defines the 
new parameter B by the transcendental equation 



B-hxB = b. 

The angular distribution function can be written then as 

oo 

W M (9,B) = ljydyJ (9y)e-y 2 / 4 exp 



(17) 







AB 1 





;is) 



The presented expansion method is to consider the term [y 2 ln(|/ 2 /4)]/4i? as 
a small parameter. Then the W M is expanded in a power series in 1/B: 



1 1 



with 



n=0 



(19) 



W n (9,t) 



1 



ydyJ,\\y\e- y21 ' 



y2 ln^)l\ (20) 



9 2 = xlB = A7in t 



Za\ 
pv J 



B(t). 



This method is valid for B > 4.5 and 9 2 < 1. 
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The first function Wo(6,t) have simple analytical form: 

^o(M = pexp(-£), (21) 



6' 



~ tint. 

t — > oo 



(22) 



For small angles, i.e. 9/9 = 9/(xcVB) = O less than about 2, the Gaussian 
( |2l| ) is the dominant term. In this region, W\(9,t) is in general less than 
Wo(9,i), so that the corrections to the Gaussian is of order of 1/B, i.e. of 



order of 10%. An alternative way of obtaining the approximate solution (21 ) 
of ([7]) for the thick targets is given in Appendix. 

A good approximate representation of the distribution at any angle is 
W o (0,t) + B~ l W\{9 } t), where 



^i(M = pexp(-g 



Ei[t)-]n(t 

9 2 \9 2 



+ 1 



(23) 



Ei(Q) = Ei(Q) + iri, and Ei(Q) = — j e is the exponential integral [15 



-e 



4 Moliere's determination of the screening 
parameters 

On the one hand, Moliere writes the elastic Born cross section for the fast 
charged particle scattering in the atomic field as follows: 



° B (X) = * R (X) 1 



F A {px) 



* R (x) q B {x). 



(24) 



For angles x small compared with a radian the exact Rutherford formula has 
a simple approximation: 

91 



° B {x) 



4irn t(l — cos x) 2 X 4 
' Q B (X). 



Q B (X) 



nn tx 



(25) 
(26) 



Here, F A is the atomic form factor and q B {x) is the ratio of actual to the 
Rutherford scattering cross sections in the Born approximation. 
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Then the screening angle Xa m the Born approximation one can represent 
via F A or q B (x) by the equations 



lim 



lim 



1 - 



F A (PX)\ 2 dx , 1 



H In q 

X 2 



Q B (x)dx 



X 



In q 



(27) 
(28) 



with an angle q such as 



Xo < ? < Uv ~ Xc 



(29) 



and the angle xo ~ m e aZ 1/ ' 3 /p. 

Moliere approximation for the Thomas — Fermi form factor F T _ F (q) with 
momentum transfer q can be written as 



F T - F {q) M = E 



CiXf 



(30) 



in which 



d = 0.35, c 2 = 0.55, c 3 = 0.10, 

Ai = 0.30A, A 2 = 4Ai, A 3 = 5A 2 . 

In the case where the Born parameter a = 0, the equation ( |27| ) for the 
screening angle can be evaluated directly, using the facts that q(0) = and 
lim q(() = 1. Then with use of (24) and (30) can also be obtained the 

(— >oo 

following approximation for [x' a ) PQ I 



(XaY = [exp(C £ - 0.5)] = VU7A XoA, 



(31) 



where A = m e aZ 1//3 /0.885. We should note that in Refs. [H E] admitted a 
misprint, namely, the factor A = 1.0825 in Eq. (31) should be replaced by 
A = 1.065 = V / Tl3. 

On the other hand, Moliere writes the nonrelativistic Born cross section 
in the form 



° B {x) 



oo 

J pdpJ (2kp sin |) *» u {p) 



(32) 
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where the Born phase shift is given in units of h = c = 1 by 



2 C U x (r)dr 

v 



\J 'r 2 — p 2 



- [ U x (r= Vp 2 + z 2 ] dz. 



(33) 



Here, k is the wave number of the incident particle, the variable p corresponds 
to the impact parameter of the collision, and U\{r) is the screened Coulomb 
potential of the target atom 



a 



Ux(r) = ±Z-A(Xr) (34) 
r 

with Moliere's fit to the Thomas — Fermi screening function A(Ar) 

A(Ar) ~ 0.1e" 6Ar + 0.55e- 1 - 2Ar + 0.35e-°- 3Ar . (35) 



In order to obtain a result valid for large a (i.e. not restricted to the Born 
approximation) and also for large angles x ( U P to 90° and beyond), Moliere 
uses in his calculations of the screening angle a WKB technique. 

We point out that the value of Xa will depend both on screening properties 
of atom, and on <r(x)-approximation used for its calculation. 

The exact formulas for the WKB differential cross section o(x) an d cor- 
responding q(x) in terms of an integral are given in Moliere's paper [1] as 
follows: 



oo 

J pdp J (fcx/>)|l - exp 



(36) 



q(x) 



(kxY 

A a 2 



oo 

J pdpj (kxp)^l - exp 



(37) 



with the phase shift given by 



kJr) — k 



dz. 



(3f 



where k r (r) is the relativistic wavenumber for the particle at a distance r 
from the nucleus, and the quantity p is seen to be impact parameter of the 
trajectory or 'ray'. As before, k is the initial or asymptotic value of the 
wavenumber. 
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When k r (r) is expanded as a series of powers of U\(r)/k, then the first- 
degree term yields the same expression for $ M (p) as Eq. (33). The Born 
approximation for (36) is obtained by expanding the exponential in (36) to 
first order in parameter a (13). 

Relations (26) and (28) between the quantities cr s (x), q B (x)i an d Xa 
remain valid also for a(x), ?(x)> an d Xa- 

Despite the fact that the formulas (36), (37) are exact, evaluation of these 
quantities was carried out by Moliere only approximately. 



To estimate (|37|), Moliere used the first-order Born shift (|33|) with (|34|) 
what 

q(x) 



and (35), what is good only to terms of first order in a, and found 

Aia(l — ia) 2 



1 



(x/xof 



1 1 , X 

— + — + lg — 

— ia ha 2xo 



0.81 + 2.21 

2 



(39) 



Here, ip is the so-called 'digamma function', i.e. the logarithmic derivative 
of the T-function ip(x) = dha.T(x)/dx. 



He has fitted a simple formula to the function 3? [^(io)] in Eq. (39) 



U [ip{ia)} 



a 4 + — + 0.13 



(40) 



Using (|40|) in Eq. (|39|) and expanding with neglect of higher orders in a 2 and 

7.2 • 10- 4 ( X /xo) 4 



(x/Xo) 2 , he got 



q(x) 



1 - 



8.85 

(x/Xo) 



1 + 2.303 a z lg 



(41) 



a 4 + a 2 /3 + 0.13)_ ' 

Moliere has calculated q(x) for different values of a and as a result has devised 
an interpolation scheme, based on a linear relation between (x/Xo) 2 an d a 2 
for fixed q: 

(X/Xof « A + a 2 B a 



Calculating the screening angle defined by 



In (Xa) 



1 

— h hm 

2 ?^-oo 



q(x)dx 

X 



1 

2 



Inxo 



dq In 



X_ 
Xo 



(42) 



(43) 



and assuming a linear relation between x\ an d a 2 , Moliere writes finally the 
following interpolating formula: 

Xa ~ XoVl.lS + SJea 2 . (44) 

Critical remarks to his conclusion see in the review of V.T. Scott [Ml . 
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5 An alternative determining the screening 
parameters 

To obtain an exact correction to the first-order Born screening angle (x' a ) B > 
we will carry out our analytical calculation in terms of the function Q e i{i]) 

oo 

Q( V ) = 2n [ a( X )[l-Jo(xv)\xdx= [ d 2 p \l - cos[A<Z>(p, ff )] , j (45) 



where the phase shift is determined by the equation 

A$(p,r/) = $(p+)-$(p_), p± = p±rj/2p. 

Substituting the expression for the cross section 



Tin t x 4 



q(x) 



into ( 45 ) , we rewrite it in the form: 



(46) 



(47) 



n t Q( V ) = 2x1 J [I- Jo(xv)\ q{x)X^d X . (48) 
o 

For the important values of rj of order 1/xc or less, it is possible to split the 



last integral at the angle q (29) into two integrals 



!(v) = j [i - Mxn)] q(x)x 3(i x 

o 



[i - Mxn)] q(x)x d x+ / [i - Mxn)} i(x)x d x 



(49) 



For the part from to q, we may write 1 — Jq(xv) = X 2? ? 2 /4, and the integral 
Ji reduces to a universal one, independently of n: 



iikv) = -j I v(x)dx/x- 



(50) 
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For the part from q to infinity, the quantity q(x) can be replaced by unity 
and the integral I2 can be integrated by parts. This leads to the following 
result for I 2 : 



77 



1 - ln(^) + In 2 -C B + 0{<;ri) 



(51) 



Integrating (50) with the account (43), substituting obtained solutions back 



into (48), and also using definition 



ln(Xc/Xa) 2 + l-2C B = ln( Xc / X :) 2 
clS cL result for Q{rf) we obtain: 

(XcV) 2 



Q(v) 



2n n t 



In 



xW 
4 



In 



Xc 
Xa 



(52) 



(XcV? ln [ V 2 (Xaf 



2n t 



4 



Finally, considering definition of 6 C (15), we can represent Q{rf) by the fol- 
lowing expression: 



(53) 



Then, the screening angle \' a can be determined via Q{rj) by a linear equation: 



ln M 2 = ln (T) 



2711]'- 



Za 

Jp 



-1 



Let us present the quantity Q e i{v) m the form: 

Q( v ) = Q B ( V )-A GO [Q( V )]. 



(54) 



(55) 



Then, on the one hand, using (53), the difference A cc [Q(rj)] between the 
Born approximate Q B {i]) and exact in the Born parameter values of the 
quantity Q (77 ) can be reduced to a difference in the quantities ln (x' a ) and 



A cc [Q(r / )]^g fl (r / )-g(r / 



Za 



In M-ln (x'a)' 
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On the other hand, this difference can be reduced to a difference Aq(x) 
Q B (x) -q(x) : 



A CC [Q( V )} = 2irJ X d X Aa{ X )[l - Mxv)} = ^ J ^Ag(x)[l - J Oct)],(56) 
o o 



and for (37) with (46), as a result of integration, (56) becomes 



A C c[Q(v)} = W 



47T17' 



(Za\ 


2 






(Za\ 


2 




I 



*ii +< f)j 



+ c 



E ( , 



where 



K [-0 (1 + ia)] = 3? (1 — ia)] = 3? (ia)] = 3? [ip (-id)] 

OO j 

= - G - + a2 E^T^) = - G . + /W. 

71=1 ' 



(57) 
(58) 



(59) 



-oo < a < oo, 

,2m- 1 



V'(l) = — C B , and f(a) = a 2 Yl™=i [ n ( n2 + a 2 )] is an 'universal function of 

a = Za/ (3 '. 

As consequence, we can get the following rigorous relations between the 
quantities In ( X ' a ) and ln {x'aY ■ 



Acc[ln W)l 



K[^(l + ia)-^(l)] , 
a 2 Y°° \n{n 2 + a 2 )}- 1 . 



(60) 
(61) 



We point out that the relations (|58j), (|60j), and (61) are independent on 
the form of electron distribution in atom and are valid for any atomic model. 

From ( |58) ) also follows an expression for the correction to the exponent 
of ([7]). Since ln[p(r/)] = —n tQ, we get: 



A cc [\ng( V )] = Hg(rj)] -ln[<f ( V )] 

2 



(62) 



47r?7 n t i 



Za 
Jp 
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For the specified value of rf = 1/xl, using the definition of Xc (15), we can 
evaluate this correction: 



Aoo [In g{Xc)\ 



xl 



Xc 



An not 



f(a) = /(a). 



(63) 



The formulas for the so-called 'Coulomb correction' (CC), defined as a 
difference between the exact and the Born approximate results, are known 
as the Bethe — Bloch formulas for the ionization losses [TH] and the formulas 
for the Bethe — Heitler cross section of bremsstrahlung |17Q The similar 
expression was found for the total cross section of the Coulomb interaction of 
compact hadronic atoms with ordinary target atoms [19] . Were also obtained 
CC to the cross sections of the pair production in nuclear collisions [201 IZE] 
as well as to the solutions of the Dirac and Klein-Gordon equations [22J. 
Specificity of the expressions received in the present work is that they define 
CC not to the cross sections of the processes, but to screening angle and an 
exponential part of the distribution function. 



6 Relative Coulomb corrections to the Born 
approximation 



Let us write ( |61[ ) in the following way 

(Xa) = {x'aY exp 



At) 



(64) 



The exact (Coulomb) relative corrections to the Born screening angle (x' a )' 
can be then represented as 



Sec (X'a) 



X'a-ix'aY _ A(xi) 



[X'a]' 



(X'a)' 



exp 



1. 



(65) 



The relative CC to the exponent g B (n) at n 2 = 1/Xc> as follows from (63), is 
also determined by this quantity, 5 cc (Xa) = $cc \d(Xc)]- Moreover, because 



AW( X c,t) = W M -W 1 



J (9r])Ag(xc)v d V, 



(66) 



accounting J °° dnn Ja(6n) = 0, we get 

1 The more complicate formal expression for CC was derived by I. 0verb0 in |18) . 



13 



Sec [W M ( X c,t)} 



AW( Xc ,t) Ag( Xc 



W B ( Xc ,t) M g B ( Xc ) 



exp 



./ 



(t) 



I. 



(67) 



Thus, 



5 CC = 5 cc (x' a ) = 5 CC [g(xc)] = 5 CC [W M (x c ,t)) 



Note also that from the Moliere equation (|39|) and to it equivalent 
8.85 



?(x) 



1 



(x/Xo) 



l + 4a z 



X 



f 



Za 



0.543 



(6* 



2 X oJ J V P 

ensues an approximate expression for the relative correction of exact to the 
Rutherford cross section 5(a) = (a — a R )/a R : 



5(a) 



8.85 

(x/xoY 



l + Aa 1 



x 

2Xc 



-/ 



Za 

T 



0.543 



(69) 



The expression in the inner brackets is close in its form to the inner parts of 
the Bethe — Bloch [16], Bethe-Maximon [T7] formulas, and also formula for 
the total cross section from [19]. 

To compare the second-order relative corrections to the first-order results 

(2) (2) A k A It 

(5 CC and 5 M ), which correspond to Eqs. (64) and (12), respectively, and 



to investigate their Z-dependence, we first present these equations in the 
following approximate form: 



(X'a) - {X'aY 



(70) 



W) - (Xa) 

The last expression follows from 



Xa*X B a \l 1 + 3.34^ «xf 



(71) 



,(72) 



where Xa = (x£)/1.080 and xl = (x«) A-080. Then, (|70f and Q becomes 



A(X^) 



1.204 



1.670 



Za 

T 



+ o 



o 



ZaX 



Za 



(73) 
(74) 
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Additionally, in order to assess the accuracy of the Moliere theory in 
determining the screening angle x' a i we a ^ so define the absolute and relative 
differences between the values of 8$ (xL) an d $cc(Xa) as we U as between the 
8^ (x^) and 8 CC (x!x) by the relations 



8 {2) (y') 

CCM\Xa) 
ScCM (Xa) 



a(2) 



(X'a) 



f(2) 



A, 



(X'a) 



8 



(2) 



r(2) _ r(2) 

u Af y CC 

8 {2) - 5 

y Af U CC 



1 - 



8, 



(2) 

CC 



CC 



(75) 
(76) 



Table 1 presents the Z dependence of the second-order corrections to the 
first-order results. It shows that the values of the relative corrections 8^1 for 
heavy target atoms with Z ~ 80 does reach 40%. From it is also obvious 
that with the rise in the nuclear charge the absolute inaccuracy Acl M {Xa) 
of the Moliere 's theory in determining the screening angle x' a increases to 
approximately 16%, and the corresponding relative error 8^cm{x'o) does not 
depend on Z and is about 28%. 



Table 1. The Z dependence of the second-order corrections defined by Eqs. 



(73)-(75) 



M 


Z 


Za 


r(2) 
°cc 


8 {2) 


A (2) 


r(2) 

U CCM 


Be 


4 


0.029 


0.001 


0.001 


0.000 


0.286 


Al 


13 


0.094 


0.011 


0.015 


0.004 


0.280 


Ti 


22 


0.160 


0.031 


0.043 


0.012 


0.280 


Ni 


28 


0.204 


0.050 


0.070 


0.020 


0.286 


Mo 


42 


0.307 


0.113 


0.157 


0.044 


0.280 


Sn 


50 


0.365 


0.160 


0.222 


0.062 


0.279 


Ta 


73 


0.533 


0.342 


0.474 


0.132 


0.278 


Pt 


78 


0.569 


0.390 


0.541 


0.150 


0.279 


Au 


79 


0.577 


0.400 


0.554 


0.154 


0.278 


Pb 


82 


0.598 


0.431 


0.598 


0.166 


0.279 



We can also calculate the exact ('Coulomb') corrections according to the 
formulas (61), (65| and compare them with the results of the calculation of 
Moliere, using (76) (Table 2). For this purpose, we must first calculate the 
values of the function f(a) = 3t[ip(l + ia)] + C E . 
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Table 2. The Z dependence of the exact corrections, defined by Eqs. 



(61), (65), and (76) 



M 


z 


Zct 


S 


f(Za) 


U CC 


^■CCM 


r 

®CCM 


Be 


4 


0.029 


1 


2012 


0.0010 


0.001 


0.000 


0.286 


Al 


13 


0.094 


1 


1928 


0.0107 


0.011 


0.004 


0.280 


Ti 


22 


0.160 


1 


1758 


0.0306 


0.031 


0.012 


0.280 


Ni 


28 


0.204 


1 


1602 


0.0487 


0.050 


0.020 


0.287 


Mo 


42 


0.307 


1 


1127 


0.1046 


0.110 


0.047 


0.297 


Sn 


50 


0.365 


1 


0799 


0.1436 


0.154 


0.068 


0.306 


Ta 


73 


0.533 





9710 


0.2758 


0.318 


0.157 


0.330 


Pt 


78 


0.569 





9467 


0.3067 


0.359 


0.182 


0.336 


Au 


79 


0.577 





9414 


0.3125 


0.367 


0.187 


0.337 


Pb 


82 


0.598 





9263 


0.3316 


0.393 


0.205 


0.342 



From the digamma series 

^(1 + a) = 1 - C E - —— + jVl)" [C(n - 1 

1 + a L 



,n-l 



\a\ < 1, (77) 



ra=2 



where ( is the Riemann zeta function, leads the corresponding power series 
for ft[V>(l + ia)] = U[ip(ia)]: 

1 oo 

K[^(m)] =1-C E -— - + J2(- 1 ) n+1 [C^n + l)]a 2n , \a\ < 2, (78) 

a n=l 

and the function f(a) = a 2 J2^=A n ( n2 + ° 2 )] 1 can be represented as 



.. oo 

/(«) = i- IT ^+E(- 1 ) n+1 [c( 2n + 1 

n=l 



,2n 



|a| < 2, (79) 



= 1 + 0.2021 a 2 - 0.0369 a 4 + 0.0083 a 6 - . . . 

1 + a z 

A equivalent way to estimate f(a) to four decimal signs, is to present the 

(80) 



sum from ( 59 ) in the following form [T7j : 

oo 

E = (i+« 2 )" + E(- a2 )""^+i)-i], 

n=l 

= (l + a 2 ) + 0.20206 - 0.0369a 2 + 0.0083a 4 - 0.002a 



Eq. (80) is sufficient to evaluate this sum up to a < 2/3 = 0.667. 
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The calculation results for (79) and (80), relative CC and their difference 
with the Moliere corrections are given in Table 2. Some results from Tables 
1,2 illustrates Figure 1. 



60 




Figure 1: The dependence of the relative Coulomb (5 C c) and Moliere 
(5^ ) corrections, and their difference (A CCM ) from nuclear charge Z. 



It will be seen from Table 2 that for the light elements up to Z = 28, all the 
exact corrections coincide with the corresponding second-order corrections 
of Table 1. Beginning with Z = 42, the relative CC are lower than above- 
mentioned 5cc (x'a) > $cc (Xa) > an d this discrepancy increases approximately 
to 10% for the heavy elements with Z ~ 80. The magnitude of A cc [in (%„)] 
from (61) is reached for Z ~ 80 the value of order of 30%. The size of 
the corresponding relative CC for these values of Z is approximately 40%. 
The absolute and relative differences with Moliere corrections at Z ~ 80 are 
about 20% and 34%, respectively. 
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Thus, in the case of scattering on targets with large Z the such cor- 
rections as 5 cc {x'a) and A CCM [Xa) become significant and should be taken 
into account in determining the lifetimes of relativistic elementary atoms in 
experiments with nuclear targets (see, e.q., [3]). 

The application of formulas (60 ), (61 ) to the Migdal theory of the Landau — 
Pomeranchuk — Migdal effect, allowed to agree the results of theory and on 
experiment, was considered in [23 J. 



7 Conclusions 

1. We have obtained the rigorous relations for the first-order Born and 
exact values of the quantities Q(rj), ln[g(r))], and x' a i which do not 
depend on the shape of the electron density distribution in the atom 
and are valid for any atomic model. 

2. We also calculated the Coulomb corrections A cc = A cc [in (Xa)] — 
&cc [Infl'(Xc)] and relative corrections 5 CC = 5 C c{x' a ) = &cc [g{Xc)\ = 
8cc \W M (xc t)] with nuclear charge ranged from Z = 2 to Z = 82 and 
showed that for Z = 28 these correction are approximately 5 percent, 
and their maximum values for Z ~ 80 comprise the order of 30 percent 
and 40 percent, correspondingly. 

3. Additionally, we evaluated absolute and relative accuracies of the Molie- 
re theory in determining the screening angle x' a - m particular, we con- 
cluded that for a target material, which intended to be used in the 
DIRAC experiment, these corrections are about 2 and 29 percents in 
the case of a Ni target (Z = 28), and they amounts to 16 and 33 
percents, respectively, for a tantalum target [Z = 73). 
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Appendix: Approximate solution for the thick 
targets 



We can also obtain the approximate solution (21) of Eq. ([7]) for the thick 
targets in the following simple way. If the total number of collisions is 



iVo = 2nn t / o(x)xdx > 1, 



then at small angles like XoV C 1, we can write 



i - Mxv) 



which will make it possible to reduce the integral Q to a much simpler one: 



W M (9,t) = j T}dr}J Q {e7})exp 
o 

In view of the that 



-2vm t— / a(x)x d X 



lim a(x)x 3 -> 0, 



the corresponding integrand from (83) is a convergent integral 



J3) 



J4) 



o{x)X Z dx < oo, 



55) 



and the integration of ( 83 ) , using the fact 



dr\r\ Jo(9rj) = 2c 



r(i) 
r(o) 



56) 



where T(x) = (x — 1)! is the Euler Gamma function [21], yields the final 
result: 



with 



W u (6,t)M^exp (-t 



P = 2irn Q t / <r(x)x 3 dx. 



(87) 
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For the Rutherford law 



at large angles, when a R (x) 3> Oq = Xo, the quantity (88) take on a value 



«» = W/^)X'* = 0=, (90) 



and the approximate solution (87), as well as (21), is not applicable. 
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